Finite nuclear size correction to the bound-electron g factor in a hydrogenlike atom 
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The finite nuclear size correction to the bound-electron g factor in hydrogenlike atoms is in- 
vestigated in the range Z—l-20. An analytical formula for this correction which includes the non- 
relativistic and dominant relativistic contributions is derived. In the case of the Is state, the results 
obtained by this formula are compared with previous non-relativistic analytical and relativistic nu- 
merical calculations. 
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I. INTRODUCTION 

Recent experiments on measuring the bound-electron g factor in hydrogenlike carbon (C 5+ ) reached an accuracy 
of about 2 • 1CT 9 The same accuracy is expected to be achieved soon for some other \ow-Z ions. To obtain the 

corresponding precision in the theoretical predictions for the bound-electron g factor, the relativistic, QED, nuclear 
recoil and nuclear size corrections must be evaluated (see |^-||] and references therein). In the present paper we derive 
a relativistic formula for the finite nuclear size correction to the bound-electron g factor in the case of an arbitrary state 
of a hydrogenlike atom. This formula provides a sufficiently accurate evaluation of the correction under consideration 
in the range Z = 1-20. In the case of the Is state, the results obtained by this formula are compared with the results 
obtained by the related non-relativistic formula jjj and with the results of the relativistic numerical evaluation jl| . 

The relativistic units (K = c = 1) and the Heaviside charge unit (a = e < 0) are used in the paper. 

II. BASIC FORMULAS 

We consider a hydrogenlike atom, placed in a homogeneous magnetic field, A(r) = [B x r\/2. We assume that the 
interaction of the electron with the magnetic field is much smaller than the fine structure splitting and much larger 
than the hyperfine structure splitting, if the nucleus has a nonzero spin. The energy shift of a level a to the first order 
in the magnetic field is 

AE a = (*„|Vb|*„), (2.1) 

where 

V B = ^(B-[?xa]). (2.2) 

Assuming that the vector B is directed along the z axis, the energy shift reads 

AE a = fiogBrrij , (2.3) 

where = ^ is the Bohr magneton, m,j is the z-projection of the angular momentum, and g is the bound-electron 
g factor. In the case of a point-charge nucleus a simple calculation, based on the Dirac equation, yields Q 

9D = , - , (k^ ^ . (2.4) 

Here n = (— iy +l+ ^(j + h) is the relativistic angular quantum number, j is the total angular momentum of the 
electron, I = j ± 1/2 defines the parity of the state, E nK is the energy of the state 

E nK = — m , (2.5) 



n r = n — \k\ is the radial quantum number, n is the principal quantum number, 7 = \J n 2 — (aZ) 2 , and 
\/K+7) 2 + (aZ) 2 . 
The finite nuclear size induces a deviation of the g factor from its Dirac value, 



1 



g = <?d + Ag. 

To find the nuclear size correction by perturbation theory, we have to evaluate the expression 



AE a = 2 



(*«|Vk|* n )(*n|Vn|*«} 



(2.6) 



(2.7) 



where Vf = V — Vc defines a deviation of the potential from the pure Coulomb one, Vc = —aZ/r. After a simple 
integration over the angular variables, the Ag value reads 



2nm \r-^\ (nn\VF\n' K)(n' K\ra x \nK) 
9 = 3 (J + 1) \f ~ 



where |n/t) is a two-component vector defined by 

|n/t) 



rg nK ,(r) 
rf nK (r) 



(2.8) 



(2.9) 



g nK and f nK are the upper and lower radial components of the Dirac wave function defined as in [Bj , a x is the Pauli 
matrix acting in the space of the two-component vectors, and the scalar product of the two-component vectors is 
defined by 



(a\b) = dr r 2 (g a g b + f a f b ) . 



(2.10) 



The sum (|2.8 )can be evaluated analytically using the method of the generalized virial relations for the Dirac equation 



proposed in 



10f . Equations (3.26)-(3.29) of Ref. |L0| can be written as 

(E n > K - E nli )(n' n\r s \nn) = ~s(n' ^iayr 3 ^ 1 \nn) 



(2.11) 



(E n > K — E nK ) (n' K,\i<j y r s \tik) = 2m(n' ' n\a x r s \nn) + s(n' ' n\r s \nn) — 

— 2k(tl n\o z r s ~ x \tik) , 

(E n i K i + E nK )(n' n\<7 z r s \nn) — 2m(n' n\r s \nn) + s{n' n\(j x r s ~ l \nn) — 

—2aZ(n , K\tj z r s ~ 1 \nK) , 

[E n i K t + E nK )(n! K\<j x r s \nn) = — s(n' K,\(j z r s ~ l \nn) + 2K(n'K|r s_1 \nn) 

— 2olZ{ti! n\a x r s ~ x \tik), 

where a x ,a y ,a z are the Pauli matrices. From these equations we obtain 



{n'n\ra x \nn) = ( k—^- - J— ) (nf k\tik) + (E nK - E n > K ) 
m z Ira 

K 



x (n k 



nr 



m \ 



E nK ) ria y + mra x + aZia y — k<t 2 

2k J 



From equation ( 2.15 ), taking into account that 



and 



(n k\tik) = S n 



\n' k) (ti'k\ = I — \nn) (nn\ , 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 
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wc find 



n'jtn , , , , i i , 

E\n K)(n k \ra r \nn) ,„ , , „, 
p IF = ( p - (™\P\™)) l««>. 



where 



^=4 



/ 771 \ 

V " K ~~ 2k/ r * (T? ' + mr(Tx + _ K(J z 



Finally, we should evaluate the expression 



2kto 



((nK\VpP\nK) — (nK\Vp\nK) (nK\P\nn)) 



(2.18) 



(2.19) 



(2.20) 



We assume that the nuclear charge distribution is described by a spherically symmetric density p(r) — p(r), which is 
normalized by the equation 



J dr p(r) 



1 . 



The Poisson equation gives 



AVp(r) = A-naZ\p{f) - S(r)}, 



(2.21) 



(2.22) 



where A is the Laplace operator. When integrated with Vf, the radial functions g(r) and f(r) can be approximated 
by the lowest order term of the expansion in powers of r. It follows we have to evaluate the integral 



/ 



1= dr rV 7_2 V P . 



Using the identity 



rP = 



(/? + 2)(/? + 3) 



Ar /3+2 



(2.23) 



(2.24) 



and integrating by parts, we find 

oo oo 

I = [ drr 2 — Ar 2 '' V F = I dr 

J 2 7 (2 7 +l) J 



1 



AnaZ 



27(27- 



OO 

— J dr r 2 r 21 p(r) 



27(27 + 1) 

aZ 



27(27+ 1) 



r 2 ^ AV F 



where 



(r 2 ^) = J dfr 2 ^p{r). 



(2.25) 



(2.26) 



For the correction to the g factor we obtain 

k 2 T(2 7 + 1 + n r )2 2 ''- 1 



j(j + 1) 7(27 + l)r 2 (2 7 + l)n r \(N - k)N 2 -t 



+2 



(n 2 + (N - k) 2 ) - 2k^) - 2n r (N - k) - 2k 



x(aZ) 2 ^(r^). 



(2.27) 



For ns-states, which are of particular interest, the expansion of this expression to two lowest orders in aZ yields 
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i + (azy {- + 



1 12n 2 -n-9 



An 2 {n + 1) 



-2*(3) - *(2 + n) 



where ^(x) = 4- logr(x). For the Is state, we have 



{r 2 \og{2aZmr / n)) 



Ag=°-(aZfm 2 (r 2 ) 



l + {aZ) 2 [ 2-C- 



: \og(2aZmr)) \ 

W) ) 



where (7=0.57721566490 is the Euler constant. In the non-relativistic limit, we find 



(2.28) 



(2.29) 



for ns states and 



A.^-(«2)V{, 2 ) 



A 9 = ^( Q Z)V{r 2 > 



for npi states. In the case of the Is state, the expression (2.3C) coincides with the related formula in |7| 



(2.30) 



(2.31) 



III. NUMERICAL RESULTS 



In Tabic 1 we compare the Ag values for the Is state obtained by formula ( 2.29| ) with the non-relativistic results 
of Rcf. (it corresponds to equation (|2.30 ) of the present paper) and with the relativistic numerical results of Ref. 
J|. To calculate (r 2 logr) in equation ( [2.29 ), we considered the homogeneously charged sphere model for the nuclear 
charge distribution. As one can see from the table, the relativistic contribution to Ag becomes comparable with 
the current experimental accuracy for ions with Z > 12. It will be also important for lower Z ions, provided the 
experimental accuracy is improved by an order of magnitude. 
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TABLE I. The finite nuclear size correction Ag for the Is state. 



z 


( r 2\l/2 

fm 


Ref. H (= Eq. (2.30)) 
[10" 9 1 

L J 


Eq. (2.29) 
[10" 9 1 

L J 


Ref. il 
[10" 9 1 

L J 


1 


0.862 


0.00003768 


0.00003770 


< 0.01 


2 


1.671 


0.0022655 


0.0022705 


< 0.01 


4 


2.390 


0.074154 


0.074741 


0.09 


6 


2.468 


0.40031 


0.40710 


0.42 


8 


2.693 


1.5064 


1.5499 


1.56 


10 


3.006 


4.5822 


4.7810 


4.78 


12 


3.057 


9.827 


10.426 


10.40 


14 


3.123 


19.00 


20.54 


20.47 


16 


3.263 


35.38 


39.05 


38.90 


18 


3.427 


62.52 


70.53 


70.28 


20 


3.478 


98.15 


113.4 


113.15 
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